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We investigate the spin relaxation induced by acoustic phonons in the presence of spin-orbit 
interactions in single electron Si/SiGe lateral coupled quantum dots. The relaxation rates are 
computed numerically in single and double quantum dots, in in-plane and perpendicular magnetic 
fields. The deformation potential of acoustic phonons is taken into account for both transverse 
and longitudinal polarizations and their contributions to the total relaxation rate are discussed 
with respect to the dilatation and shear potential constants. We find that in single dots the spin 
relaxation rate scales approximately with the seventh power of the magnetic field, in line with a 
recent experiment. In double dots the relaxation rate is much more sensitive to the dot spectrum 
structure, as it is often dominated by a spin hot spot. The anisotropy of the spin-orbit interactions 
gives rise to easy passages, special directions of the magnetic field for which the relaxation is strongly 
suppressed. Quantitatively, the spin relaxation rates in Si are typically 2 orders of magnitude smaller 
than in GaAs due to the absence of the piezoelectric phonon potential and generally weaker spin- 
orbit interactions. 

PACS numbers: 03.67.Lx, 71.70.Ej, 72.25.Rb, 73. 21. La, 73.22.Dj, 85.35.Gv 



I. INTRODUCTION 

For more than a decade semiconductor quantum dots 
have been in research focus for quantum information 
processing^— The experimental control over the elec- 
tron spin in quantum dots has seen enormous progress, 
with lateral gated GaAs structures demonstrating the 
state of the art^ Magnetic^ and electrio 9 - coherent spin 
manipulations have been demonstrated, while the elec- 
tron spin measurement has been achieved using spin 
to charge conversion techniques. 10 For the latter, dou- 
ble dots 1 - 1 - have proven especially useful exploiting the 
spin Coulomb blockaded However, the GaAs spin qubit 
seems to have reached its fundamental limit for the co- 
herence which is due to the nuclear spins inherent in all 
III-V semiconductors J£r— 

Materials composed of atoms without nuclear magnetic 
moment, such as Si and C, seem a natural solution for the 
problem of the nuclear induced decoherence That is 
why Si based quantum dots have recently seen a revived 
interest. Although the quantum dot technology is not yet 
as mature as in GaAs, several perspective setups are be- 
ing actively pursued ,22r— We note that a spin to charge 
conversion was reported recently.— In addition to the ab- 
sence of nuclear spins, Si seems potentially advantageous 
because of weaker spin-orbit interactions, promising less 
decoherence, and a stronger g factor, allowing spin con- 
trol at smaller magnetic fields. 

On the other hand, the electron effective mass in Si 
is larger than in GaAs, so Si dots must be smaller at a 
given orbital energy scale. In addition, and perhaps more 
seriously, a major issue for silicon based quantum com- 
putation is the valley degeneracy of its conduction band 
electrons i 28 ' 29 In the bulk, the conduction band minima 
are located at the X valleys, that is at k v ~ 0.84fco, 
v = 1, ...6, toward the six X points of the Brillouin 



zone, where k$ = 27r/ao and ao = 5.4 A is the lattice 
constant^ In a (OOl)-grown Si heterostructure the valley 
degeneracy is partially lifted due to the presence of the 
interface and/or due to strain^ leaving a twofold con- 
duction band minimum, the ±z valleys, which are sepa- 
rated from the fourfold excited valley states by at least 
10meV ) 29 ' 30 i 32 large enough to neglect the upper four 
valleysi 28 ! 29 ' 33 The remaining twofold valley degeneracy 
is lifted if the perpendicular confinement is asymmetric. 
Then the orbital wave functions become symmetric and 
antisymmetric combinations of the single valley states^ 
which are separated by the energy difference called the 
ground-state gap^ 3 - (or valley splitting). 

In recent years the origin and possible control of the 
valley splitting has been in focus. Measurements in sili- 
con heterostructures reveal a valley splitting of the order 
of /xeV i 31 i 35 ~ — On the other hand, theoretical estimates 
of perfectly flat structures propose a splitting about three 
orders of magnitude larger; 41 Taking into account de- 
tailed properties of the interface (e.g., roughness), exper- 
iment and theory come to an agreement ) 33 ' 34 > 42 ~ — and 
additional (in-plane) confinement allows the valley split- 
ting to reach values of the order of rneV^ In Si/Si02 sys- 
tems, the splitting can even be tens of meV^ 9 - - — Here we 
assume that the splitting is at least 1 meV and we can use 
the effective single valley approximation ; 29 ' 52 in which 
only the lowest valley eigenstate is considered. This 
choice is strengthened by the fact that electron spins in 
valley-degenerate dots would not be viable qubitsi 28 ' 29 ' 52 

In the single valley approximation, the Si dot resem- 
bles the fairly well understood GaAs one. The main goal 
of this article is to carry out a comparison on a quantita- 
tive level, providing realistic values for the electron spin 
relaxation as available for GaAs dotsi^r— 

The relevant sources of spin relaxation in GaAs quan- 
tum dots in a magnetic field are the electron-phonon cou- 
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plings modeled by the piezoelectric and deformation po- 
tential theoryi 56 ' 58 Since silicon is not piezoelectric, only 
the deformation potential mechanism remains. From this 
point of view, a Si dot is closer to InGaAs than to GaAs, 
as in InGaAs the relative importance of the deformation 
versus the piezoelectric potential is enhanced due to a 
larger g factor^ However, differently from (In) GaAs, in 
which transverse acoustic phonons do not contribute to 
the deformation potential coupling, both transverse and 
longitudinal acoustic phonons cause spin relaxation in 
Sii££r— j n .j. ne p resen t WO rk special attention is given to 
the importance of the transverse phonons and to the role 
of the dilatation and shear potential constants, S<j and 
S u respectively, which parameterize the electron-phonon 
coupling strengths. 

In Si, the spin relaxation/decoherence rates were com- 
puted perturbatively for single dot single electron 6 ^ and 
single 66 and double 67 dot singlet-triplet transitions. Also 
the spin relaxation due to the modulation of electron g 
factor by the phonon-induced strain was investigated. 68 
Experimentally the rates were measured on quantum dot 
ensembles 69,70 and on a many-electron quantum dot ) 71 ' 72 
and a few electron quantum doti^l We obtain the spin 
relaxation rates non-perturbatively, using exact numer- 
ical diagonalization, for a wide range of magnetic fields 
and interdot coupling. Our results allow us to discuss 
regimes beyond the validity of perturbative treatments 
and to specify the accuracy of common approximations. 

The article is organized as follows. In Sec. II we de- 
fine our model of the double quantum dot, the electron- 
phonon interaction, and the spin relaxation. In Sec. Ill 
we briefly review the Si single and double dot spectra, 
paying attention to the states' orbital symmetries. In 
Sec. IV. A we investigate the spin relaxation in a single 
dot in an in-plane magnetic field. Comparing with a re- 
cent experiment, we find that (i) the results of the experi- 
ment indicate that the main spin relaxation channel is the 
mechanism we study here, and (ii) the spin-orbit coupling 
strength of ~ O.lmeVA seems realistic for the Si/SiGe 
lateral quantum dots. We also present analytical formu- 
las for the spin relaxation rate in the lowest order of the 
spin-orbit interactions. Comparing with exact numerics, 
we demonstrate that these formulas are quantitatively 
reliable up to modest magnetic fields of 1-2 T. We find 
that a further analytical simplification often adopted, the 
isotropic averaging of the interaction strengths, leads to 
a result correct only within an order of magnitude. In 
Sec. IV. B we deal with the double dot case and demon- 
strate that here the spin relaxation rate is sensitive to 
the spectral anticrossings (spin hot spots) ^ especially 
if the magnetic field is perpendicular to the heterostruc- 
ture. For in-plane fields, the anisotropy of the spin-orbit 
interactions leads to the appearance of easy passages — 
magnetic field directions in which the relaxation rate is 
quenched by many orders of magnitude ! 56 ' 58 These re- 
sults are analogous to the GaAs quantum dots. Finally, 
we conclude in Sec. V. 




FIG. 1: Orientation of the double dot in the coordinate system 
x — [100], y = [010]. The potential minima, sketched by two 
circles, are parameterized by the position vectors ±d or by 
the distance 2d and the angle 8. The in-plane magnetic field 
orientation is given by the angle 7. 

II. MODEL 

Within the two-dimensional, effective mass and ef- 
fective single valley approximations a laterally coupled, 
top-gated double quantum dot (DQD) in a silicon het- 
erostructure with the growth direction z = [001] is de- 
scribed by the Hamiltonian 

H= ^- + V(x,y) + ^ B a-B + H hl + H d , (1) 
2m 2 

with m the effective electron mass and r = (x, y) is the in- 
plane position vector. The vector of the kinetic momen- 
tum KK. = — iftV + eA, with e the proton charge, consists 
of the canonical momentum and the vector potential in 
symmetric gauge, A = B z (—y/2,x/2). We neglect the 
orbital effects of the in-plane magnetic field component, 
a good approximation up to roughly 10 T for the usual 
heterostructures.— The in-plane coordinates are chosen 
to be the crystallographic axes, that is, x = [100] and 
y = [010]. We use the biquadratic mode l 29 ' 75 for the 
electrostatic confinement potential, 

V(x, y) = X -mul min{(r - d) 2 , (r + d) 2 }, (2) 

where ±d denote the positions of two potential minima. 
Alternatively, the minima are parameterized by the in- 
terdot distance 2d and the angle between the main dot 
axis d and x, denoted by S. Characteristic scales are 

given by the confinement energy E$ = fkuo and the con- 

1/2 

finement length Iq — (h/mujo) . The Zeeman energy 
in Eq. JT]) is proportional to the effective Lande factor g 
and the Bohr magneton /ib, while cr — (a Xl a y ,a z ) is the 
vector of the Pauli matrices. The magnetic field is given 
by B = (_B|| cos 7, £?|| sin 7, B z ), where 7 is the angle be- 
tween the in-plane component of B and x. The geometry 
is plotted in Fig. Q] 

Extrinsic spin-orbit coupling (SOC) leads to additional 
terms in the Hamiltonian of two-dimensional systems 
without inversion symmetry* 5 - Structure inversion asym- 
metry arising, for example, due to an electric field ap- 
plied along the growth direction, results in the Bychkov- 
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Rashba Hamiltonian 



Hbr = a (a x K y - cr y K x ) , 



(3) 



with an electrically tunable coupling parameter a. The 
bulk inversion asymmetry of zinc blende semiconductors 
such as GaAs, is not present in silicon with diamond 
structure and the corresponding linear and cubic Dres- 
selhaus interactions are absent. However, a silicon het- 
erostructure is of C2 V symmetry, so there is also a gener- 
alized Dresselhaus term t 76 ' 77 



H d =f3 {-<7 X K X + GyKy) , 



(4) 



which is identical to the Dresselhaus Hamiltonian stem- 
ming from the bulk inversion asymmetry in III-V semi- 
conductors. The parameter j3 depends on the interface 
(step height) and well width (number of Si-atomic lay- 
ers). An alternative parameterization utilizes the spin- 
orbit coupling lengths ?br = ft 2 /2ma and l d = h 2 /2mf3. 
In this work we assume an asymmetric quantum well 
along z with a and /3 of comparable strengths^ 

Our parameters for the numerical calculations are 
as follows: The system of interest is a (OOl)-grown 
SiGe/Si/SiGe quantum well, where the thin Si layer is 
sandwiched by the relaxed SiGe with a Germanium con- 
centration of about 25%. The bulk electron effective mass 
of the X valleys is anisotropic in the direction transverse 
and longitudinal to the corresponding fc„-vector, given 
by mi — 0.916m e and nit = 0.191m e , respectively, where 
m e is the free electron mass. The in-plane mass of the 
z valley states is therefore the transverse mass. Due to 
the tensile strain (in-plane) of the Si layer the effective 
mass is slightly increased^ 8 - compared to the unstrained 
bulk Si and we use m = 0.198m ef 2& The effective Lande 
factor is g = 2 42 ' 71 and the SOC strengths are set to be 
a = 0.05 me VA and f3 = 0.15 me VA, respectively^ Our 
choice is based on results of theoretic tight-binding cal- 
culations of Ref. [77], as experimentally the SOC in silicon 
dots has not been measured up to date. We use the con- 
finement energy Eq = 1 meV equivalent to the confine- 
ment length Iq = 20 nm, which corresponds to realistic 
dot sizes i 71 1 79 

We consider transitions mediated by phonons, with the 
accompanying spin- flip allowed due to the presence of the 
spin-orbit interactions. In a (OOl)-grown quantum well, 
the electron-phonon coupling for intravalley scattering is 
the deformation potential of transverse acoustic (TA) and 
longitudinal acoustic (LA) phonons given by£2 r 64 ' 80 




2pVc x Q 



bkxj QR -bQ,xe- iQR ], (5) 



with 



(6) 



where Q = (q, Q z ) is the phonon wave vector, Q is 
its unit vector and R = (r, z) is the electron posi- 
tion vector. The summation includes all polarizations 81 



A = TA1, TA2, LA, and c\ is the corresponding sound 
velocity. The phonon annihilation (creation) operator is 
denoted by b Qy) and the polarization unit vector reads 
e. The mass density is given by p and V is the volume 
of the crystal. The deformation potential strength is set 
by the dilatation and shear potential constants Sd and 
S u , respectively. Note that unlike in GaAs there is no 
piezoelectric phonon potential in Si. 

We define the single-electron spin relaxation rate, 
which is the inverse of the spin lifetime T\, as the sum 
of the transition rates from the upper Zeeman split 
ground state, called Tg in the following, to all lower-lying 
states 'J j. with opposite spin. Each individual transition 
rate is evaluated using Fermi's Golden Rule in the zero- 
temperature limit^ i 56 ' 75 



where = (Fg|e l( 3' R | , I'^) is the matrix element for the 
corresponding initial and final states and fiw-1-4. is the en- 
ergy difference between these states. We evaluate Eq. ([7]) 
numerically using the parameters p — 2.3 x 10 3 kg/m , 
and Ct = 5 x 10 3 m/s for TA, and q = 9.15 x 10 3 m/s for 
LA phonons. The choice of deformation potential con- 
stants is not unique^ 2 - We use 3d = 5 eV and 3 U = 9 eV 
according to Ref. HH, noting that other combinations 
such as (3 d ,H u ) = (1.1, 6.8) eV& (1.13, 9.16) eV, and 
(— 11.7,9)eV— appear alternatively. The needed elec- 
tron wave functions and energies are obtained numeri- 
cally as eigensystem of the Hamiltonian in Eq. (TIJ , which 
we diagonalize with the method of finite differences us- 
ing Dirichlet boundary conditions. The magnetic field is 
included by the Peierls phase and the diagonalization is 
carried out by the Lanczos algorithm. Typically we use a 
grid of 50x50 points, which results in a relative precision 
in energy of 10 -5 in zero magnetic field. 



III. SINGLE ELECTRON STATES 

A. Energy spectrum in zero magnetic field 

In order to understand the details of the spin relaxation 
in silicon DQDs, we review briefly their electronic prop- 
erties in zero magnetic field, including group theoretical 
classification, the influence of SOC and the most impor- 
tant quantities for experiments. The Hamiltonian Eq. (fTJ) 
for B = and without SOC has C2 V ® SU(2) symmetry. 
We can then label the orbital states according to the ir- 
reducible representations T^, i — 1,...,4 of the Abelian 
point group C2 V noting that each state is doubly spin- 
degenerate due to SU(2). This is done in Fig. [2j where 
the energy spectrum vs. the interdot distance in units of 
lo is plotted. Note that the potential V, Eq. ([2]), was cho- 
sen such that the states converge to Fock-Darwin states^ 
in the limit of zero or infinite interdot distance. In the 
following we focus on the intermediate region where the 
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FIG. 2: (color online) Calculated energy spectrum of the Si 
double quantum dot with respect to the interdot distance at 
zero magnetic field. The states are labeled (colored) according 
to the irreducible representations Ti of C2v On the right- 
hand side we give the highest orbital momentum of associated 
single dot states (Fock-Darwin states). The tunneling energy 
T is also shown. 

interdot distance is comparable to the confining length. 
This is typically the region of experimental interest, as 
well as the one in which numerics becomes indispensable. 
Here we find several level crossings which may be lifted 
in the presence of SOC. Such anticrossings, also called 
spin hot spots ^ are of great importance for spin relax- 
ation as we will see later. However, the linear SOC terms 
[Eq. ([3]) and Eq. (jlj] do not lead to level repulsion in the 
first order although allowed by symmetry^ We conclude 
that in zero magnetic field the DQD spectrum of silicon 
does not exhibit relevant spin hot spots. 

For many applications including quantum dot spin 
qubits, the important physics happens at the bottom of 
the spectrum. We denote the spin-degenerate ground 
state as T\ = T$ and the first excited state as T2 = Ta to 
indicate the symmetry under inversion I xy . The energy 
difference between these states is parameterized by the 
tunneling energy^ T = (E& — E$) /2, a characteristic 
quantity for DQDs directly measurable experimentally. 86 
Note that within the single valley approximation we as- 
sumed a valley splitting of at least 1 meV which exceeds 
2T at all interdot distances. 

Using a linear combination of single dot orbitals 85 
(LCSDO) we can approximate the exact wave functions 
by analytical expressions. Let \& n) /(r) be the Fock- 
Darwin state (omitting spin), where n is the principal 
and / the orbital quantum number- 8 - Then the lowest or- 
bital eigenstates of the DQD can be approximated using 
the Fock-Darwin states centered at the potential minima 
as 

r s = W + [* 0l o(r + d) + *o,o(r-d)], 

T A = iy_[* ,o(r + d)-*o,o(r-d)]. (8) 

Here N± are normalization constants. Calculating the 
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FIG. 3: Calculated tunneling energy vs. interdot distance for 
zero magnetic field calculated by exact numerical diagonal- 
ization (dotted line), exact LCSDO formulas (thin solid line) 
and leading order approximation (Eq. ©, thick solid line). 
The tunneling energy for a finite perpendicular magnetic field 
(B z = 2 T, dashed line) is given for comparison. 

eigenenergies as the expectation values of the Hamilto- 
nian, Eq. ([T]), for zero magnetic field and without SOC, 
we obtain the tunneling energy as plotted in Fig. [3] It is 
in excellent agreement with the exact numerical result. 
In the limit of large interdot distances the leading order 
reads 

T^E ^-e-W^\ (9) 

which is a good approximation if 2d/lo > 2.5. 

In principle SOC terms affect the tunneling energy. 
However, it was showrt 8 ^ that this correction is of fourth 
order in the spin-orbit strengths a and/or (3. For our 
parameters here it is of the order of peV and therefore 
negligible for all experimental purposes. 



B. Energy spectrum in non-zero magnetic field 

In a perpendicular magnetic field without SOC the 
group of the Hamiltonian becomes the Abelian point 
group C2 . The only remaining symmetry operator is the 
total inversion I xy and the one-dimensional irreducible 
representations have either symmetric or antisymmetric 
base functions. The spectrum of a DQD in the perpen- 
dicular magnetic field is plotted in Fig. [4] The Zeeman 
interaction lifts the spin degeneracy and the ground state, 
denoted as Eg, is spin-polarized. Up to a certain mag- 
nitude of B z (about 1.5 T for 2d/l = 2.5; see Fig. gj, 
the first excited state is r|, and the spin relaxation is 
the transition between these two wave functions with the 
same orbital parts and opposite spins. For larger mag- 
netic fields the Zeeman splitting exceeds the orbital exci- 
tation energy and the first excited state is P^, which has 
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FIG. 4: Calculated energy spectrum of the DQD with inter- 
dot distance 2d/lo = 2.5 plotted against the perpendicular 
magnetic field. The thick line indicates Ig, the lowest state 
with opposite spin-polarization as the ground state. 



the same spin polarization as the ground state. For even 
higher fields more states fall below r|, which all con- 
tribute to the spin relaxation. Note that the level spac- 
ings of interest at moderate magnetic fields are smaller 
than the assumed valley splitting which again justifies 
the single valley approximation. 

Within the LCSDO the single dot wave functions ac- 
quire a phase when shifted. The building blocks in 
Eq. © now read * n ,j(r ± d) exp[±iei3 z r • (z X d) /2h] 
and we can repeat the computation of the tunneling 
energy with the result plotted in Fig. [3] One can see 
that the perpendicular magnetic field reduces the tun- 
neling energy. This can be understood qualitatively as 
the renormalization of the confinement length, which is 
replaced by the effective (magnetoelectric) confinement 
length Is, where = Iq^ + I^ 4 " with the auxiliary quan- 

1/2 

tity l 6 = (2H/B z e) 1 . The tunneling energy simplifies in 
the limit of large interdot distances to 



T^E - 



-d 2 (2i" 2 -i 2 B /- 4 ) 



(10) 



Note that for B z = we have Ib = lo and we obtain the 
results of Sect. IIII Al 



IV. SPIN RELAXATION 

A. Single Quantum Dot 

Before we proceed to DQDs, we first discuss a sin- 
gle quantum dot in an in-plane magnetic field, B = 
Bu (coS7, sin7, 0), which already features anisotropics 
and relaxation rate spikes due to spin hot spots as we 
will see. Removing the linear SOC terms in Eq. (|T|) by a 
unitary transformation ; 8 ^ 88 the relaxation proceeds due 
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FIG. 5: (color online) Spin relaxation rates of a single electron 
in a single quantum dot vs. in-plane field for 7 = 135°. The 
total rate (solid black line) and its contributions of the TA 
phonons (solid red line) and the LA phonons (solid blue line) 
are shown. The spin hot spot at By ~ 8.3 T causes the spin 
relaxation rate to increase up to the orbital relaxation rate 
(dash-dotted line). The three dots give the experimental data 
of Ref. |7l| fitted by a B 7 curve (dotted line). 



to a spin-orbit-induced effective magnetic field 



56.58 



B 



off. 



-B\\ 



cos 7 sin 7 



^br 



Id 



( sin 7 

V ^br 



cos 7 



(ii) 

which is perpendicular to the external magnetic field. 
The matrix element Mf-j. in Eq. (0 is proportional to this 
effective magnetic field, which results in the spin relax- 
ation rate being proportional to the squared and inverse 
effective spin-orbit coupling length L , 56 ' 58 i 89 



L~ 



7-2 , 7-2 



2sin(27)/ b - r 1 /J 1 , 



(12) 



which is anisotropic since it depends on 7. However, the 
anisotropy disappears if one of the SOC lengths is dom- 
inant, particularly for /3 = 0. Thus, experimental ver- 
ification of anisotropic single dot spin relaxation would 
verify the existence of the generalized Dresselhaus term, 
Eq. (j4|). The anisotropy is strongest if Zbr = Id, with 
the maximal rate at 7 = 135° and the minimal rate at 
7 = 45°. 

Figure [5] displays the numerical results for the spin 
and orbital relaxation rates with 7 = 135°. The TA 
and LA phonon contributions to the total spin relaxation 
rate (solid black line) are given to clarify their relative 
importance. We find that the rate for magnetic fields 
up to about 5 T essentially results from the TA phonons. 
They do not depend on Sd since the scalar product in 
L>q, Eq. ©, vanishes. An important observation is the 
strong enhancement of the total spin relaxation rate at 
S11 « 8.3 T. This is due to a spin hot spot which appears 
at a point at which the Zeeman splitting is equal to the 
level spacing of the Fock-Darwin states. The anticrossing 
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induces a strong mixing of the spin states which abets 
spin relaxation. The spikes appear with equal height for 
any in-plane field orientation 7 as here the rate is given 
by the orbital relaxation rate^ (dash-dotted line), which 
is independent of 7. 

In Fig. [5] we also plot the spin relaxation rate as mea- 
sured in Ref. [Til- First, the observed power dependence 
corresponds to the coaction of spin-orbit interactions and 
deformation phonons The energy conservation forbids 
a direct electron-nuclear spin flip-flop in finite magnetic 
fields. This process becomes allowed if accompanied by 
the emission of a phonon, yielding a relaxation rate pro- 
portional to B 5 £Q Second, the order of magnitude agree- 
ment indicates that our choice of the spin-orbit strength 
is realistic, even though a direct fitting is not possible 
(the angle 7 was not reported and the dot was not in a 
single electron regime). 

We now derive analytic formulas for the spin relax- 
ation rate valid for weak in-plane magnetic fields. Treat- 
ing the spin-orbit coupling perturbatively, we are able 
to evaluate Eq. ([7]) analytically. The total rate, r sp i n = 
is given by the contributions (A' = TA, LA) 



^TA 
spin 



spin 



spin 



m 2 Zg 
A ' 2A-K P cl,h w 




(13) 
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FIG. 6: Ratio of the relaxation rate of the approximations 
and the exact numerics vs. in-plane magnetic field. The low 
B- field limit, represented by Eqs. (|14p and (115[l . is shown 
by the solid lines, and the numerically evaluated isotropic 
average approximation, Eq. (|16|l . by the dashed lines. The 
contributions of the TA and LA phonons are given in the top 
and bottom panel, respectively. The magnetic field is in-plane 
and 7 = 135°. The constant line at 1 (dotted) is a guide to 
the eye. 



The energy parameter T>\, reads 



for the transverse, and 



^TA — — >-h, 



35 



and 



T) 2 — ~ 2 X n n 
^LA — + 5" d " u 



3 „ 
35" 



(14) 



(15) 



for the transverse and longitudinal branches, respectively. 
The weak versus strong magnetic field limit is determined 
by the conditions £ \ <C 1 and £ \ 3> 1 respectively, where 
£\ = gnBBlB/(fic\)£& Here, the crossover £ \ = 1 is 
found at B» = 1.4 T for transverse, and B\\ = 2.6 T 
for longitudinal acoustic phonons. Comparing with the 
exact numerics, we find that the error of the value of 
Eq. (H) is less than 10% up to B\\ = 0.8 T for TA, and 
up to B\\ = 2 T for LA phonons. In any case, the error is 
less than 5% if By < 0.5 T. 

The integral in Eq. can be done analytically only 
exceptionally, such as in the single dot case. Therefore, 
one often employs isotropically averaged deformation po- 
tentials to simplify the treatment! 62 ' 63 ' 91 This amounts to 
average Dq, Eq. ([6]), over phonon directions distributed 
uniformly in three dimensions, 



D, 



Ql 



Here it leads to Eq. (05} with 



D Q \ 2 dn. 



^TA.iso 



2«> : 



_~ 2 



15 



(16) 



(17) 



LA,: 



W) 2 



—2 



(18) 



for the longitudinal contribution to the total rate. For 
our choice of parameters, we get (L>q A ) = 3.29 eV and 
(Z)L A ) = 8.44 eV. Comparing Eqs. (17]) and JTS) with 
Eqs. (fL4]l and (flip)) , we find that the averaging, Eq. (fTBT) . 
leads to rates which are 2.3 (TA) and, for our parameters, 
1.4 (LA) times larger than the actual rates. Note that 
if we use instead of Eq. (|T6|) the averaged deformation 
potentials as reported in Ref. l63l . we obtain relaxation 
rates that differ in the low-i?-field limit from Eqs. (fLT| 
and ((15|) by a factor of 3.4 for the TA, and 3.2 for the 
LA contribution. 

The isotropic average approximation becomes exact if 
the matrix element \M^\ is independent of the phonon 
direction. However, this directional invariance is not ful- 
filled in lateral dots which are strongly anisotropic in the 
perpendicular versus the in-plane direction. To assess 
the quality of the approximations, we compare the corre- 
sponding relaxation rates with the exact numerical result 
in Fig. |6] for magnetic fields up to 10 T. Our measure is 
the ratio between the rate of the approximation and of 
the numerics, which we plot for the TA and LA contribu- 
tions separately. The parameters in Fig. [5] are identical to 
Fig.[SJ We find that the analytical results (solid lines) de- 
viate significantly from numerics for fields beyond the low 
B-field limit. The especially large discrepancy at around 
8.3 T, where the ratio is close to zero, stems from the fact 
that the analytic approximations assume no level cross- 
ings of the initial state. Thus, it accounts for neither spin 
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FIG. 7: (color online) Calculated spin relaxation rate in a 
DQD with tunable interdot distance in a perpendicular mag- 
netic field. The rate is given in inverse seconds by the color 
with the scale on the right. The j/-axis is calibrated in interdot 
distance (left) and tunneling energy at B = (right). 



FIG. 8: (color online) Calculated spin relaxation rate in a 
DQD as a function of the interdot distance and the orientation 
of the in-plane magnetic field (By = 4T). The rate is given 
in inverse seconds by color with the scale on the right. The 
y-axis is calibrated in interdot distance (left) and tunneling 
energy at B = (right). The dots main axis is along [100]. 



hot spots, nor the transition into excited states. Equa- 
tions (fl"7]) and (fl8|) result in curves parallel to Eqs. (fT4|) 
and (jTSJ) , but shifted by the discrepancy factors 2.3 and 
1.4 for the TA and LA contributions (not shown). Nu- 
merical evaluation of the spin relaxation rates via Eq. ([7]) 
using the average of Eq. (fT6|) leads to a discrepancy rep- 
resented by the dashed line. We find that even in highly 
anisotropic (2D) lateral dots, the discrepancy factor is 
only of the order of 1. It is therefore expected to be 
legitimate to use the isotropic averaging also for more 
complicated dot geometries, such as the double dot, or 
a biased dot, where it can lead to significant simplifica- 
tions. 



B. Double Quantum Dot 

We now move to a DQD case, where we take 2d as a 
variable parameter, noting that it could stand for cither 
the actual separation between two dots, or a gate-tunable 
coupling between dots of fixed distance. From the experi- 
mental point of view, it is more convenient to characterize 
a double quantum dot via the tunneling energy. We plot 
our results with respect to the interdot distance and give 
also the corresponding tunneling energy at zero magnetic 
field computed numerically (Fig. [31 dotted line). 

The spin relaxation rate in a DQD as a function of both 
the interdot distance in units of Iq and the perpendicular 
magnetic field is shown in Fig. [7] We find that the plotted 
area is dominated by the spikes which come from spin hot 
spots and that there are no easy passages, that is there 
is no possibility for a fixed magnetic field to change the 
interdot distance from zero to infinite without passing 
through any of these peaks. For small fields, here B z < 
3 T, we have only one relevant spike, which comes from 
the anticrossing of Tg and T A (see Fig. @|. For larger 



magnetic fields, crossings with higher orbital states occur 
which may but need not lead to spin hot spots, depending 
on the symmetry of the crossing states. 

If the field is applied in the plane, the spin relaxation 
depends on the orientation of the magnetic field with 
respect to the crystallographic axes because of the inter- 
play between Bychkov-Rashba and generalized Dressel- 
haus SOC. Once we rotate the coordinate system by S 
around z, the effective Zeeman field, Eq. (|11|) . read o 56 ' 58 

Bf =-5|i [l£ cos( 7 - 6) - 1 sin (7 + S)] + 

+ B {] y [l£ sin( 7 -5)- if cos( 7 + S)] , (19) 

where the tilted axes are such that x is parallel to d. 
Since the first excited orbital state Ta transforms like 
x, only the first term in Eq. (|19j) can lead to spin hot 
spots for moderate magnetic fields in the intermediate 
regime. As an example, we plot in Fig. [8] the spin re- 
laxation rate in a DQD aligned along [100] (S = 0°) in 
a magnetic field Bu = 4 T varying the orientation 7 and 
the interdot distance. A sharp peak occurs at 2d/lo = 1.8 
but is intermittent at 7c = 18° where the first term of 
Eq. (fT9|) vanishes. Note that this angle, defining the easy 
passage, depends on the SOC lengths and is thus sample 
and setup dependent. An experimental determination of 
the easy passage angle would provide information about 
the relative SOC strengths via the relation tan7 c = a//3. 

Equation (fT9| shows that the easy passage depends 
also on the DQD orientation with respect to the crys- 
tallographic axes. We note that for a DQD with the 
main axis along the [110] direction (6 = 45°), the corre- 
sponding field orientation of the easy passage is universal 
(7c = 135°), independent of the SOC strengths. We plot 
the spin relaxation rates for this case in Fig. |9] for com- 
pleteness. 
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FIG. 9: (color online) Calculated spin relaxation rate in a 
DQD as a function of the interdot distance and the orientation 
of the in-plane magnetic field (By — 4T). The rate is given 
in inverse seconds by color with the scale on the right. The 
y-axis is calibrated in interdot distance (left) and tunneling 
energy at B — (right). The dots main axis is along [110]. 

V. CONCLUSIONS 

We investigated the spin relaxation of a single electron 
confined in a quantum dot in a laterally gated Si/SiGe 
heterostructure. We considered the spin relaxation to be 
an inelastic transition in which the spin-flip is enabled 
by the presence of the spin-orbit interactions, while the 
energy difference between the initial and final state, aris- 
ing from the applied magnetic field, is taken away by 
an acoustic phonon. We studied relaxation rates varying 
the interdot coupling from strong (a single dot regime) 
to negligible (a double dot regime), as well as varying 
the strength and orientation (in-plane or perpendicular) 
of the magnetic field. 

We adopted the single valley, effective mass, and two 
dimensional approximations, within which our results are 
numerically exact. Whereas the latter two are known 
to be well justified for lateral quantum dots, the sin- 
gle valley approximation breaks down once the valley 
splitting drops below the orbital energy scale, ~meV, 
and additional states appear in the lowest part of the 
dot spectrum. Concerning the spin relaxation, however, 
these states are irrelevant as the matrix elements for the 
phonon induced intervalley transitions are greatly sup- 
pressed. This is so because long wavelength phonons, 
which arise due to the small transition energy (q < 0.1 
mil -1 at 3 T), are ineffective in coupling states with dis- 



parate Bloch wavefunctions (k v w 10 nm 1 ); see Eq. (2) 
in Ref. and the discussion therein. 

We found that the spin relaxation in Si dots is roughly 
comparable to that in GaAs, although it bears certain 
differences. Namely, in the single dot the relaxation rate 
in Si is proportional to B 7 , being due to the deformation 
phonon potential, in contrast to the B 5 dependence in 
piezoelectric GaAs. We compared our theory with ex- 
perimental data, which confirm the magnetic field power 
dependence and show that the spin-orbit strengths of the 
order of O.lmeVA are to be expected in Si/SiGe quan- 
tum dots. We also derived an analytical expression for 
the relaxation rate treating the spin-orbit interactions in 
the lowest order. We find it an excellent approximation 
to the numerics up to magnetic fields of 1-2 T. A further 
simplification, the isotropic averaging, makes the analyt- 
ical result to differ from the exact one by a factor of the 
order of 1. 

We showed that in the double dot the relaxation rate 
is a much more complicated function of the magnetic 
field and the interdot coupling, the two parameters most 
directly controllable experimentally. This is due to the 
fact that the rate is strongly influenced by spin hot spots, 
which occur at much lower magnetic fields in the double 
dot compared to the single dot. The anisotropy of the 
spin-orbit interactions leads to the rates dependent on 
the magnetic field direction with respect to the crystal- 
lographic axes. In a double dot, where the rotational 
symmetry of the potential is broken but the reflection 
symmetry is preserved, this anisotropy results in the ap- 
pearance of easy passages — special directions of the ex- 
ternal magnetic field which assure a strong suppression 
of the relaxation rate. From these directions the ratio 
of the spin-orbit strengths can be found. Compared to 
GaAs, in Si the easy passage position relates directly to 
the linear spin-orbit strengths without being influenced 
by the spin-orbit interaction cubic-in-momenta. 

Finally, we observed that compared to GaAs, the spin 
relaxation rates in Si are typically 1-2 orders of magni- 
tude smaller, as a result of the absence of the piezoelec- 
tric phonon interaction and generally weaker spin-orbit 
interactions. 
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